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^ . Possible saturation of betatron acceleration of dust particles behind strong shock fronts 

pL^ ' from supernovae is considered. It is argued that the efficiency of the nonthermal dust 

destruction should be substantially lower than the value estimated from a traditional 
description of betatron acceleration of dust grains behind radiative shock waves. The 
inhibition of the nonthermal destruction can be connected with the mirror instability 
developed in the dust component behind strong shocks with the velocity 3 times exceeding 
. the Alfven speed. The instability develops on characteristic time scales much shorter the 

(-H ' age of a supernova remnant, thus its influence on the efficiency of dust destruction can be 

Q^. substantial: in the range of shock velocities 100 km s^^ < < 300 km s^^ the destruction 

^ ' efficiency can be an order of magnitude lower that normally estimated. 

^ ; I. INTRODUCTION 

At present it is commonly accepted that dust grains are efficiently destroyed behind 
\ shock fronts from supernovae [1, 2]. A rough estimate of the dust destruction rate in the 

\ Galaxy can be obtained from the following simple consideration [2]: the total mass of the 

CN ■ dust destroyed in unite time is = vsNiJ^^/Ci where vsN is the Galactic supernova rate, 

r] is the efficiency of dust destruction by a shock with a given velocity, M is the mass of 
supernova shell with a given velocity; the kinetic energy of the shell is assumed by order 
(<~^ ' of magnitude to be equal to the explosion energy E ~ Mv^ /2. It is commonly assumed 

that dust grains most efficiently are destroyed by the shock waves with the velocities 
V > 100 km s^^ — at lower velocities the destruction efficacy drops exponentially [1]. The 
fraction of destroyed dust behind the shock with the velocity fs,o.4 = 100 km s~^ is around 
rj ~ 0.4, hence the total rate of the dust mass decrease in the Galaxy is estimated as 



X 



Md ~ 0.8usN^§^ (1) 

C<0.4 

where Esn = 10^^ erg is the explosion energy, ^ = 100—200 is the mean gas-to-dust mass 
ratio. The characteristic time of dust destruction is 



Td,d = = 1-25^^ ~ 10° year, (2) 

where Mg = 5 x 10^ Mq, the gas mass in the Galaxy, the suparnova rate is assumed to be 
i^SN = 1/100 year~^. 

It is obvious that this value of Td^d is unacceptably short in comparison with the time 
of dust formation (replenishment) in the Galaxy. Indeed, the total rate of dust injection 
into the interstellar medium (ISM) from stellar sources (including supernovae) is <3 x 
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10~^^ Mq pc~^ year""^ [3-5]. Assuming now the total dust mass in the Galaxy ~ 5x lO''' Mq, 
one can estimate the replenishment rate is t^j ^ 10^ years. In such conditions the total 
dust mass in the Galaxy must be an order of magnitude lower than the observed value, i.e. 
~ 5 X 10^ Mq. 

Moreover, it is well known that the abundances of refractory elements in gas phase of 
the ISM are always substantially lower than solar values. In other words such elements are 
"depleted" in gas phase. This depletion is connected with their condensation on the surface 
of dust grains. Even in the diffuse rarefied phase (n ~ 0.05 cm~^) the depletion factor 
can be as high as 90% for Fe, Ca, Ti, i.e. the fraction of these elements in gas phase is 
10% or lower, while in more dense regions (n ~ 3 cm~^) the fraction of Fe in gas phase is 
around 1%, and the fraction of Ca and Ti is onlu 0.1% [6,7]. This means that destruction of 
dust (and as a consequence removal of metals from condensed phase into gaseous) should 
be considerably slower process than their formation, i.e. condensation of metals on to the 
solid (dust) phase. In particular, for the diffuse ISM the following interrelation between 
the characteristic dust destruction t^^cL and formation t^j times: r^j = O.lr,^^^. In general, 
the condition t^j <^ T^^d has to be always fulfilled, though at present the mechanisms able 
to maintain this condition arc unknown. 

One of the possibilities is connected with that a considerable dust mass is confined in 
dense atomic and molecular clouds. When a shock wave falls on to a dense cloud its velocity 
inside the cloud decreases as Vs,c ~ Vs,i\/ pi/ pc, where indices c and i belong to the cloud 
and the intercloud medium, respectively. Even for HI clouds with density only one or one 
and a half orders higher than in the intercloud medium the shock velocity decreases by 
factor 3-5 , so that dust remains undestroyed behind such shock waves. Therefore, only 
dust in the diffuse intercloud gas is destroyed by supernovae shocks. With this circumstance 
accounted the above estimate can be rewritten as follows 



Md ~ Q.SiXiUsN^^. (3) 

S^s,0.4 

Here /ij = Pi{l — fc)/ Pcfc is the mass fraction contained in the intercloud medium, /c is 
the volume filling factor of clouds: 



/. = iTV-/^2^dM, (4) 

where the cloud mass spectrum dN/dM oc M^^'^ in the mass range 10"^ Mq < M < 
10^ Mq is assumed [8, 9]; the estimate gives pi ~ 0.1. As a result, the dust destruction 
time increases as ~ 2x10^/^1 and becomes comparable to the Galactic dust production 
time. At the same time the necessary interrelation Td, f <C Td^d remains unexplained. 

Apparently, one has to assume that in such conditions condensation of metals on to 
solid (dust) phase should occur directly in the ISM. A low abundance of deuterium in gas 
phase of the local interstellar medium seems to favor this conclusion. It is argued in [10] 
that the observed relative abundances of deuterium and oxygen indicate that deuterium is 
frozen-out on dust grains: the characteristic time for deuterium to freeze out is estimated 
in [10] to be around several Myr in the conditions of diffuse interstellar clouds. 

The latter means apparently that dust grains escape destruction by supernovae shocks. 
Indeed, all known data about the depletion of deuterium belongs to the nearest vicinity of 
Sun restricted by distances of 100-150 pc. The characteristic interval between successful 
explosions of supernovae in this volume is around 5 x 10^ year [11] - an order of magnitude 
shorter than the time of deuterium freezing on the dust surface. Therefore deuterium can 
be depleted in gas phase of the local ISM provided only that SNe shocks do not destroy 
dust grains. 
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In this paper we present arguments favoring that plasma instabihties behind shock waves 
prevent an efficient acceleration of charged dust particles by betatron mechanism, and thus 
inhibit their destruction. An essential circumstance here is that the dust component can 
change its dynamics only on times corresponding to the inverse dust cyclotron frequency 
ujd- This requires the presence of plasma oscillations with frequences close to (jJ^ in the 
whole spectrum, which seems plausible only under description of the dust as a separate 
dynamical component. Therefore, in this paper we describe instabilities in the interstellar 
plasma behind strong shock waves caused by the presence of the charged dust component. 

In Sect. 2 wc formulate the equation of motions of a three-fluid (ions, electrons and 
charged dust) interstellar plasma, in Sect. 3 we derive the dispersion equation for small 
oscillations in such a plasma, in Sect. 4 the dispersion equation is analyzed and the 
conditions for instabilities are discussed. Sect. 5 summarizes the results. 



An efficient destruction of dust particles behind shock fronts is determined by 
their betatron acceleration in the course of plasma compression on radiative stages [2, 
12]. Namely, during a radiative cooling and an accompanying compression of plasma 
magnetic field behind the shock front grows as i? oc p, what in turn leads 
to a proportional increase of the square of the transversal momentum of dust 
particles p\ oc B. Consequently the efficiency of dust destruction increases by 
more than an order of magnitude for the shocks with velocities 100 km s~^ < 
< Vg < 300 km s~"^. In these estimates the structure of the post-shock flow, in particular, 
the magnetic field being parallel to the front, is assumed implicitly to be stable and 
invariable over the whole post-shock region. At the same time, one can think that growing 
plasma instabilities are able to violate the symmetry of the flow and moreover to make 
the flow irregular. As a result, plasma regions with a decreased magnetic field and the 
correspondingly decreased dust velocities can arise, which result in an inhibition of the 
dust destruction efficiency. 

One of the instabilities which can cause violation of a regular magneto-hydrodynamic 
flow is the mirror instability of an anisotropic plasma with the transversal temperature 
exceeding the parallel one: T±_ 2> Ty. In the conditions behind shock fronts anisotropy of 
the dust component emerges naturally: after crossing a shock front with a parallel magnetic 
fleld a dust grain gains the transerval velocity in the restframe connected with the front 
close to the shock velocity, while the longitudinal velocity remains equal to a much lower 
velocity of random dust motions in the ISM. Below we consider a possibility for the mirror 
instability to develop in such conditions. 

Let us describe the interstellar plasma by the equations [13] 



where a = i, e, d corresponds to the ions, electrons and dust particles. Conservation of the 
momentum is described by the equations 



II. EQUATIONS OF MOTION 



dtUa + V • {nana) = 0, 



(5) 



dtinaUa) + V ■ {UaVLaUa) = 





which should be complemented by the equations 




(7) 
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and 

Vx B = ^^Z„n«v«, (8) 

as well as the equation for cur IE; here Z^. = —1, Zj = +1 and = —Zfj, so that in 
typical conditions with negative dust grains the quasi-neutrality equation has the form 
ni = ne + ZdUd- 

Let us consider plasma behind the shock in the comoving (moving with the plasma) 
restframe. Let us consider also the configuration with magnetic field poins parallel to z- 
axis (along the shock front), while the wavevector has the components k = (A;^, 0, A;||). 
In this work we neglect the effects connected with perturbations of the front. Therefore 
our results formally are valid far from the front - stability of the shock front in such a 
multi-component plasma deserves a separate consideration. 

Let us linearize the system (5)-(8) in projections for each plasma components. We 
will restrict consideration by low-frequency electromagnetic waves with the phase velocity 
satisfying the inequality 



VTi<€,Uj/k <^VTe, (9) 

where vtc and VTi are thermal velocities of the electrons and ions, respectively. For 
linearized equations equations of motion of the electrons one obtaines 

dvx e k±Tf, We 

= Ex-UJeVy-l , (10) 

mm m Ueo 

dvy e ^ ,^ , 

— = --Ey+u;eV., (11) 

= -^E^-^bEl^^ (12) 
m m Ueo 

where Ug = eBo/mc is the electron cyclotron frequency, Tg is the electron temperature. 
In equation (12) the electrons are assumed to be inertialess along the field Bq, which 
corresponds to the right hand side of the inequality (9). Then the continuity equation for 
the electrons turns out to be equivalent to the equation for a perturbation of the electron 
density derived from the latter equation: 



rieO felloe 

The ion equations of motion are written as 



dt M 



— Ex+LOiVy, 



dVy _ e 



Ey -LUiVx, 



at M 



(14) 
(15) 
(16) 
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where terms with thermal pressure are neglected in virtue of the left inequality (9). The 
continuity equation for the ions couples their density perturbations with variations of both 
longitudinal and transversal velocities: 



-(^1 + 

dt \ riio J dx 



+ 



dz 



0. 



(17) 



Contrary to the electron and ion components dust relaxes to thermal equailibrium 
much slower because of low frequency of collisions. Therefore one can assume that in the 
region where the elctrons and ions are isotropic the dust component is anisotropic. In such 
conditions in the equations of motion of the dust the anisotropic pressure tensor has to be 
explicitly used. The linearized gradient of the pressure tensor is [14] 



d 



VP 



B„ 



dp', d , 



dx 

d_ 
dx 



-DO. 



+ 



ip\\ -p±] 



+ 



dz 



(18) 



Here = ridoT^], P± = n(ioT±, nao is dust density, T]| and T± its temperature along and 
perpendicular to the magnetic field, respectively, p'j_ are small deviations of pressure 
components from the equilibrium values. It has to be mentioned that since the field line 
is determined by superposition of the external and the wave field, the pressure and p± 
in different points should have different projections on to x and z axes. Then the motion 
equations of dust get the form 



dwx 

~~dt 



— Ex - uJdWy - z 

rrid rndUdo 

ik\\ ^ ^ Br 



(19) 



dwy Ze 

= Ey + UdWx 

dt rud 



ik\\ ^ ^ Br 

rUd r>o 



(20) 



dwz 
dt 



^^^_!^(r,,-r^)^ 

rrid rud " -Bo 



ik\\ 



mdUdo 



p\l, 



(21) 



where u>d = ZeBo/nidC is the dust cyclotron frequency. 

In order to determine the perturbations of pressure pj| and p'j_ we use the following 
equations 



dp\, 

-^+P||V-w + 2p||V||w;^ = 0, (22) 
dv' 

-|^ + 2pxV- w-pxV||«;^ = 0, (23) 

equivalent to the Chew-Goldberger-Law (CGL) invariants, as well as the dust continuity 
equation 
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d f nij\ dwx dwz „ , 



dt \nda J dx dz 

It has to be mentioned that the pressure components p± in equations (22) and (23) 
are connected with both variation of density and perturbation of temperature - the latter 
having adiabatic character in CGL theory. As a result, the equation for curlB in linear 
approximation is written as 

(V X B)x,y = —{nioVx,y - neoVx,y - ZndQWx,y). (25) 

III. THE DISPERTION EQUATION 
A. Derivation of the dispersion equation 

The equations (10)-(12), (14)-(17), (19)-(25) and (7) represent a fuU set of equations 
of our problem. In general, this system of equations described the two types of low- 
frequency electromagnetic waves: Alfven and magnetosound. Our goal is to demonstrate 
that behind the shock wave mirror instability develops in the dust component. It is well 
known that mirror instability is connected with magnetosound electromagnetic field with 
the wavevector complanar to the external magnetic and the wave magnetic fields. According 
to our choice k± points along x-axis. Hence only the field components Ey, Bx (and B^) are 
connected with the magnetosound wave, while the components E.^ are By small. In such 
conditions the solution of the motion equations (lO)-(ll), (14)-(15) and (19)-(20) for the 
Fourier components ~ expi(A:xa; + k\^\^z — Lot) in the drift approximation (co <C Wi, Wd) can 
be written as 



E c ( k± 

Vx = c-^, Vy = -— Ex - —Ez ] , (26) 
Bo Bq \ fell 



= c^, = --^ (Ex + z^Ey ] , (27) 



Ey 

Wx = C— T", 

5n' 



'-(ex- ^^Ey] - ^(Ti, - T^)Bx - (28) 



V. = -i^Ez, wz = ^(7|| - T^)^ + — (29) 

here the expression for Vy is obtained with accounting (13), while in the second equation 
of (28) the term contained the field component E^ is omitted as its contribution into the 
dispersion equation is small over ZM/rud- 

The equation for curl E gives the connection between the field components of interest 

U! k\\ 
Ey = —B„Bx = -^B,. (30) 
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Prom equations (17) and (24) we have 



riio oj 

and 



(31) 



rid ^ k±wa, + k\\W;, 

UdO ^ 

With using (27), (29) and (30), it is convenient to rewrite equation (31) in the form 

In order to derive the perturbation of dust density let us write the equations following 
from (22), (23) for the Fourier components: 

V\\ = n.oT|| f 3^ - 2!^) , p'^ = UdoT^ + (34) 

II " V "do BqJ \ndo BqJ 

here equations (24), (29) and (30) are accounted. Then with accounting (28), (29) and (34) 
one obtains from (32) 



nd 



/. fen 37|| -T, \ / / k^3T 



^-Tizr" ■ (35) 



ndo Bq \ ijp- vfid ) / \ ^ ""^d 

It is obvious that T^/rrid behind the shock front has the order of the square of the 
shock velocity, while T||/md characterizes the square of the thermal velocity of the dust 
component before the shock and is therefore negligibly small. Therefore the latter equation 
can be written as 



nd B^ 



klT\ 



- o + ^ ■ (36) 
Udo Bq \^ nid^^^ ) 

The field can be determined with using the quaineutrality equation = r;^ + Z-rTd from 
(13), (33) and (36): 



.neo 



Mu? knc^^ B 



^. = -i— 2 .2 2 i^' (37) 

riio e - /cj^cl Bq 

where = (riioTe/neoM)^/'^ is the dust ion-acoustic speed. With using the latter equation 
let us write Vy as 

Ej: .UeO UJ^ k±C^gBz , , 

'"y = 9 PTo TT- (38) 

Bq riio - kicj ooi Bq 



The velocity Wy can be found from (28), (34) and (36) with accounting the above mentioned 
interrelation between the longitudinal and transverse dust temperatures: 



-Bo V Wrf 



3 + B,. (39) 



Substituting (27), (38) and (39) into the y-projection of the curlB we arrive at the 
dispersion equation 



fcV? = I 1 



Pi k\c^ 



Here V\ = Bl/A-Kp, V"^ = T^/md, p = pi + Pd, Pi = Mtim, pd = mdUdo- 

B. The mirror instability 

In the absence of dust {pd = 0) equation (40) is biquadratic 

u;^-k\Vl + c',)u;' + eqvic', = 0, (41) 

and describes fast and slow magnetosonic waves. When dust is present equation (40) is 
in general bicubic. Similar equation accounting kinetic effects is obtained and analyzed in 
details in [16]. We consider here the most simple and obvious limiting case when the wave 
propagates almost perpendicular to the magnetic field: fey <C A; « /c^ which corresponds 
to a long-wavelength modulation of the magnetosonic waves along the external magnetic 
field. Then with accounting the condition a;/A;|| S> Cg equation (40) simplifies and becomes 

uj'^-Auj^ + B = 0, (42) 

where 

A = k^ (yj + + 3^V^^ , B = (vicl + S^V^'i - 
The solution of (42) is 

u'^ = ^(^A± ^A^ - 45) . (43) 



The sign "minus" corresponds to the slow magnetosonic wave. Unstable aperiodic solutions 
for such waves exist if 5 < which corresponds to the nonequality 
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Immediately behind the shock where radiative losses are unimportant the parameter V'^, 
which characterizes the post-shock transversal dust temperature, is connected with the 
shock velocity by the relation = 9Vg/32, where Vg is the shock velocity; here we assumed 
explicitle that a dust particle crossing the shock moves with respect to plasma with the 
velocity 3us/4. The electron temperature immediately behind the front is Tg = SnieV^ /IQ, 
i.e. V'^/cl = 3mi/2me. Hence as might be thought the condition (44) fulfils always as 
mi/nie > pI Pd- In reality, however, V^^/c^ = 3/2 is valid up to the region where radiative 
cooling sets on, since betatron acceleration responsible for an efficient dust destruction 
establishes far behind the front, where the electron-ion plasma is already isothermal. 
Therefore, the condition (44) apparently does not fulfil. 

In fact, it has to be taken into account that in the domain of radiative cooling all 
parameters included into (44) vary. Namely, due to radiative cooling and an increase of 
magnetic field, and caused by it betatron acceleration the orbital velocity of dust particles 
V grows as y (X \fB., or for a frozen-in magnetic field as F oc y^. This gives for a 
cooling plasma where pT 2± const F^/c^ oc T"^; similarly in the radiation domain we have 
V|/Cg oc r~^. Therefore, in the domain of radiative cooling the inequality (44) can be 
written in the form 



where Tq is the plasma temperature immediately behind the shock. Here immediately the 

shock front we explicitly used the relation F|/c^ — 20/M^, where Ma = i's/Vao is the 
Alfven Mach number before the shock. Substituting here Tq/T > 10 and a typical for the 
ISM value p/ Pd — 100, one can find the condition for the mirror instability Ma > 3, which 
is deliberately valid for the shocks able to provide dust destruction. 



In this paper we have considered a possibility for the mirror instability to develop in the 
dust component behind the shock. Flows behind the fronts of sufficiently strong shock waves 
{Ma > 3) turn out to be unstable against the mirror instability. Consequently, a laminar 
magnetohydrodynamic flow behind the shock front supporting betatron acceleration of 
dust particles becomes violated. One can expect that it causes weakening of the betatron 
acceleration. Indeed, the mirror instability on nonlinear stages is known to result, in 
particular, in formation of magnetic voids, holes and growing loops swelling out behind the 
front to a few dust gyration radii [17], which in our conditions is r^d ^ 0.3 pc for typical 
parameters of dust particles and the shock velocity 100 km s^^. As a result, charged 
dust particles are expected to be expelled into the domains of weak magnetic field and 
decelerated as p\(x B. 

It follows from (43) that the characteristic time of the instability has an order 



where A is the perturbation wavelength. Assuming the wavelength to be equal to the 
thickness of the adiabatic SN remnant A ~ AR ~ O.lRg and p/pa equal to a typical in the 
ISM value 100, one can find to be of the order of the remnant age. This estimate is 
an upper limit of the growth time - short-wavelength perturbations grow obviously faster. 
The perturbation wavelength is restricted from below by the kinetic limit and by the size 




(45) 



IV. CONCLUSION 




(46) 
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of the nonuniformity of the post-shock flow. The latter corresponds to the characteristic 
cooling length. The radiative cooling time can be estimated as Tc ~ fcT/An ~ 10^^ s for 
temperature T ~ 3 x 10^ K, typical for the beginning of the radiative stage, and density 
n ~ 1 cm~^; the corresponding scale is A ~ 0.3 pc. The size of the remnant with the 
expansion velocity ~100 km s~^ in a medium with mean density n ~ 1 cm~^ is 30 pc, 
hence the lower estimate of the characteristic time of the mirror instability on these stages 
is of 1% the remnant age. It can be expected therefore that the instability being developed 
in the whole possible wavelength range will efficiently restrict the rate of dust destruction. 
As mentioned above that at nonlinear stages the mirror instability results in breaking up of 
a regular flow with magnetic field parallel to the front, and leads essentially to formation 
of loop-like domains with a decreased magnetic field - magnetic voids [17]. It allows to 
estimate the factor of weakening of dust destruction assuming that magnetic field behind 
the shock does not increase due to radiative cooling, but in average remains rather the 
same as if the compression due to cooling would be absent. The corresponding decrease of 
the destruction efficiency would become an order of magnitude in the range of the "danger" 
velocity interval 100 km s~^ < Vg < 300 km s""^. Firm estimates can be obtained only from 
numerical modeling of the instability. 
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